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In this paper we look at the gravitational spin-spin interaction between macroscopic 
astronomical bodies. In particular, we calculate their post-Newtonian orbital effects of 
order 0(c~^) on the trajectory of a spinning particle with proper angular momentum s 
moving in the external gravitomagnetic field generated by a central spinning mass with 
proper angular momentum J. It turns out that, at order 0{e) in the orbiter's eccentricity, 
the eccentricity the pericenter and the mean anomaly rates of the moving particle are 
affected by long-term harmonic effects. If, on one hand, they are undetectable in the 
Solar System, on the other, maybe that in an astrophysical context like that of the 
binary millisecond pulsars there will be some hopes of measuring them in the future. 

1. Introduction 

In recent years the topic of tlie general relativistic post-Newtonian effects of or- 
der 0{c^'^) on various features of the motion of spinning particles freely orbiting 
a central astronomical body has received great attention both theoretically and 
experimentally. More precisely, when a spinning particle moves in an external grav- 
itational field one has to describe both its spin precession and the influence of the 
spin on its trajectory^. 

The geodetic, or De Sitter, precession^ refers to the coupling of the static, grav- 
itoelectric part of the gravitational field due to the Schwarzschild metric generated 
by a central, non-rotating object to the spin of a particle freely orbiting around it. 
It has been measured for the Earth-Moon orbit, thought of as a giant gyroscope, in 
the gravitational field of the Sun with 1% accuracy^. It should be measured for four 
superconducting gyroscopes in the gravitational field of the Earth by the important 
GP-B mission^ at a claimed relative accuracy level of 2 x 10^^. Finally, it might be 
possible to detect it also in some binary pulsar systems'^. 

The Lense-Thirring drag of the inertial frames^ is an effect due to the station- 
ary, gravitomagnetic part of the gravitational field of a rotating central mass with 
proper angular momentum J on the geodesic path of a freely falling test particle, 
i.e. considered to be not spinning itself^. In 1998 the first evidence of the Lense- 



''Such effect could be thought of as a spin— orbit interaction between the spin of the central object 
and the orbital angular momentum of the test body. For some spin— orbit effects induced by the 
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Thirring effect in the gravitational field of the Earth has been reported* with a 
claimed accuracy of almost 20%. It is based on the analysis of the laser-ranging 
data of the LAGEOS and LAGEOS II geodetic satellites. The launch of the pro- 
posed LAGEOS-like LARES satellite^ could allow to measure such effect with an 
accuracy probably better than 1%. Another interesting gravitomagnetic effect of or- 
der 0{c~^) on the orbit of a test particle has been recently derived in reference^"; it 
is due to the temporal variability of the Earth's angular momentum. Unfortunately, 
it is too small to be detected with Satellite Laser Ranging. 

The spin of the central object affects also the spin of a particle freely orbiting 
it in a way discovered by Schiff^^ in 1959. The detection of this subtle precessional 
effect, in addition to the geodetic precession, is one of the most important goals of 
the GP-B mission*; the claimed accuracy amounts to 1%. 

In this paper we are interested in looking for some orbital effects due to the spin- 
spin gravitational coupling on the geodesic path of a spinning extended particle with 
mass m and proper angular momentum s freely orbiting around a central body of 
mass M and proper angular momentum J. The gravitational spin-spin coupling in 
the quantum mechanical domain has been treated in references'^' 

The paper is organized as follows. In Section 2 we derive the gravitational Stern- 
Gerlach force with some simplifying assumptions. In Section 3 we work out the 
long-term orbital effects of such interaction on the Keplerian orbital elements of 
the moving particle. In Section 4 we look at the Sun-Mercury system and to PSR 
B1259-63 and PSR B1913-I-16 in order to see if the predicted effect could be mea- 
sured. Section 5 is devoted to the conclusions. 



2. The gravitomagnetic Stern Gerlach force 

In the context of the linearized gravitoelectromagnetism it turns out that a gravit- 
omagnetic dipole moment for a gyroscope of spin s is'^ 

— > ^ 

and the energy of interaction with an external gravitomagnetic field Bg is 
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Since, in general, the gravitomagnetic field is position-dependent, a Stern-Gerlach 
force arises from Equation (0) 
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Let us consider a central spherical rotating body as source of the gravitomagnetic 
field, so that 
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rotation of the Sun on the orbital angular momentum of the Earth-Moon system see reference'^. 
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where G is the Newtonian gravitational constant, J is the unit vector along the 
proper angular momentum of the central body and f is the unit position vector. 
By assuming J = Jz in an inertial frame with the {x, y} plane coinciding with the 
equatorial plane of the central mass, Equation (|^) becomes 
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Equation (|^)- Equation agree with the expression 
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of Equation (22) in reference^^ and with the first term of Equation (2) in reference^. 



3. The orbital effects 



Now, let us work out the long-term orbital effects on the orbit of a spinning par- 
ticle. We will consider it as an extended spherical body^ of radius I and (slowly) 
spinning with an angular velocity a, so that s — ^mPa. We will consider the grav- 
itomagnetic spin-spin force of Equation (^- Equation (^ as a small perturbation 
of the Keplerian monopole and we will adopt the Gauss perturbing equations 
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where R, T, N are the projections of the perturbing acceleration onto the radial f , 
transverse t and out-of-plane n directions of an orthonormal frame comoving with 



''It should be noticed that, according to the authors of reference^*, the straightforward extension of 
the spin-spin gravitational interaction for elementary particles to macroscopic extended rotating 
bodies should not be justified. Then, they follow a phcnomcnological approach not based a priori 
on General Relativity. 
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the orbiter, a, e, i, Q, to and M are the orbiter's semimajor axis, eccentricity, in- 
chnation, longitude of the ascending node, argument of perigee and mean anomaly, 
respectively. Moreover, p = a(l — e^), f is the true anomaly and n = V GMa~^ is 
the Keplerian mean motion. 

Now we will derive the components i?, T, N from Equation (||)- Equation (Q). 
In order to simplify the problem we will assume to consider only equatorial orbits. 
Moreover, we will further assume that the proper angular momentums of the two 
body are aligned, i.e. s = sA: as well. This is a situation quite common, e.g., with 
the planets and their moons in the Solar System to which our attention will be 
drawn. Indeed, even if of order C(c~^), the investigated effects are very small due 
to the dependence on r"'*, as can be noticed from Equation (||)- Equation (0) and 
by considering that x, y, z are all proportional to r; so, only astronomical bodies 
could yield a chance for measuring them. 

Let us make the x axis coincide with the line of the nodes: in this case 

X — cos([x) + f)f + cos i sm[uj + /)< + sin i smiw + /)n, (15) 

y = — sin(cj + f)r + cos i cos(lj + J)i + sin i cos(a; + /)n, (16) 

z — — sin i t + cos i n, (17) 

and 

x = rcos(w + /), (18) 

y = r cosi sin(cj + /), (19) 

z = r sinz sin(a; + /). (20) 

With Equation ([l^)- Equation (|20|), Equation (||)- Equation (0) become, in the 
case of equatorial orbits 

R = -^cos[2(c^ + /)], (21) 

T - -^sin[2(c^ + /)], (22) 
TV = (23) 



where s — — — ^Pa. From Equation (23) and Equation (fll)- Equation (|l2|) it 



can be noticed that the inclination and the node, which, on the other hand, is 
not defined for equatorial orbits, are not affected by the gravitomagnetic spin-spin 
force. 

By calculating Equation (pT|)- Equation (|3|) on the Keplerian unperturbed orbit 

r-Kep = -± 24 

1 + e cos / 

neglecting the terms of order 0{e^), inserting the obtained result into Equation (||)- 
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Equation (14) and averaging over an orbital revolution^ by means of 

= ^7T ^\2^/' (25) 

n(l + ecos jY 

we obtain that only the eccentricity, the perigee and the mean anomaly experience 
long-term, harmonic perturbations with half the period of the perigee. 
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4. Is it possible to measure the gravitomagnetic spin spin force? 
4.1. A Solar System scenario 

In order to measure the small effects predicted by Equation (|2^)- Equation ( [27| ) 
the most suitable system seems to be the Sun-Mercury one. In Table 1 we quote 
the relevant orbital parameters of it. For them we have used the reference From 



Table 1. Orbital parameters for the Sun and Mercury 



Parameter 


Description 


Value 


Units 


Mq 


Mass of the Sun 


1.9891 X 10^=* 


g 


R@ 


Mean equatorial radius of the Sun 


6.96 X 10"' 


cm 


£0 


Tilt of the Sun's spin axis to the ecliptic 


7.25 


deg 


^0 


Sidereal rotational period of the Sun 


24.65 


days 


A.U. 


Astronomical Unit 


1.4959787066 x IQi^ 


cm 


Mm 


Mass of Mercury 


3.302 X 10^6 


g 


Rm 


Mean equatorial radius of Mercury 


2.4398 X 10** 


cm 


Qm 


Semimajor axis of Mercury 


0.3870 


A.U. 


era 


Eccentricity of Mercury 


0.2505 






Inclination of Mercury's orbit to the ecliptic 


7.0048 


deg 


£m 


Tilt of Mercury's spin axis to the ecliptic 


0.01 


deg 




Sidereal rotational period of Mercury 


58.646 


days 


Pm 


Orbital period of Mercury 


0.2408 


years 


G 


Newtonian constant of gravitation 


6.67259 X 10"* 


— 1 ^ _9 

g cm'^ s 


C 


Speed of light in vacuum 


2.9979 X 10^0 


cm s~^ 



Table 1 it can be inferred that, if we choose the ecliptic plane as {x, y} plane, the 
assumptions made by us are satisfied rather well for the Sun and Mercury. Indeed, 
the Sun's spin axis is only slightly tilted to the vertical with respect to the ecliptic 

■^In the averaging process we consider the vector s as constant, i.e. we neglect the De Sitter and 
Lense— Thirring precessions of the spin of the orbiting particle. 
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plane, the inclination of Mercury's orbit to it amounts to few degrees and the spin 
axis of Mercury is vertical, so that it can be considered aligned to that of the Sun, 
at least to an approximate extent. The values of Table 1 yield for the Mercury's 
perihelion shift 



duj \ 
'dt/ 



= -(7.6 X 10"^° mas/y) X sin 2w. (29) 

Mercury 

Unfortunately, this is a value which is completely undetectable. 
4.2. A pulsar scenario 

Let us try to see if the predicted effects could be measured in an astrophysical 
context. The binary systems in which at least one member is a pulsar neutron star 
seem to be the optimal choice due to the large values of the eccentricity of the 
orbits of a not negligible fraction of them (70 out of more than 1300 known, i.e. 
about 5%). Moreover, the rapid rotations of the pulsars and the relatively small 
separations between the members of such systems could be helpful. 

As can be inferred from Table 2, obtained from the dataf] in reference^^, the PSR 



B1259-63 system seems to be a good choice (See, e.g., http:/ /www .jb.man.ac.uk/ 
pulsar/ for a quick outlook of the discovered pulsars and their relevant orbital 
parameters). Of course, in that case our calculations should be considered as very 
preliminary order-of-magnitude estimates, just to understand if looking at pulsars 
more thoroughly and carefully could be fruitfully or not. In the case of PSR B1259- 

Table 2. Orbital parameters for the PSR B1259-63/SS 2883 binary system 



Parameter 


Description 


Value 


Units 




Mass of the eompanion 


10 


Mq 


Rc 


Equatorial radius of the companion 


6 


Rq 


^max 


Break-up velocity of the companion 


480 


km s"'^ 


m 


Mass of PSR B1259-63 


1.4 


Mq 


r 


Radius of PSR B1259-63 


10 


km 


T 


Period of rotation of PSR B1259-63 


47.7620537 


ms 


P 


Period of revolution of PSR B1259-63 


1,236.7238 


days 


i 


Inclination of PSR B1259-63 


36 (144) 


deg 


e 


Eccentricity of PSR B1259-63 


0.86990 




^ — a sin i 

f- 


Projected semimajor axis of PSR B1259-63 


1296.4 


s 



63 the companion is SS 2883, a lOth-magnitude star of spectral class B2e, almost 
10 times more massive than the pulsar which orbits it in 3.4 years along a very 
eccentric path. In this case we will assume as {x, y} reference plane the plane 
of sky which is a fixed plane normal to the line-of-sight. For the definition of the 

"^The value of the pulsar radius r is an estimate based on the first calculations by Oppenheimer 
and Volkoffl''. 
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various orbital parameters ol the pulsar systems see reference^^. It should be noticed 
that in applying our results to PSR B1259-63 the following caveat hold. Neither 
the orbit of the pulsar can be considered equatorial|, nor the spin J of SS 2883 is 
vertical^^. Moreover, we can calculate s for the pulsar, but we cannot say anything 
about its direction. On other hand, the value can be inferred for J is plausible 
because the stars of the same type of SS 2883 are well known. 

By using the values of Table 2 we obtain for the periastron advance 

= -(3.3 X 10"" deg/y) x sin2w. (30) 

PSR B1259 63 

Notice that the period of the secular precession of the perigee of PSR B1259-63 
amounts to 1.9x10® year^ so that over reasonable observational time spans of 
few years the predicted harmonics would be similar to secular trends. According to 
Table 1 of reference^^, the experimental sensitivity amounts to 10"® deg/y for (u), 
so that also in this case the gravitomagnetic spin-spin interaction turns out to be 
too small to be detected. 

A more favorable situation occurs for the well known binary pulsar system of 
PSR B1913+16 in which the companion of the pulsar is probably another neutron 
star. The relevant orbital parameters are in Table 3 and have been taken from 
reference^^. In this case we do not know anything about the period of rotation of 

Table 3. Orbital parameters for the PSR B1913+16 binary system 



Parameter 


Description 


Value 


Units 


mc 


Mass of the companion 


1.3873 


Mq, 


m 


Mass of PSR B1913+16 


1.4411 


Mq, 


r 


Radius of PSR B1913+16 


10 


km 


T 


Period of rotation of PSR B1913+16 


59.029997929883 


ms 


P 


Period of revolution of PSR B1913+16 


3.22997462736x10-1 


days 


i 


Inclination of PSR B1913+16 


47 (133) 


deg 


e 


Eccentricity of PSR B1913+16 


0.6171309 




^ a sin i 


Projected scmimajor axis of PSR B1913+16 


2.341759 


s 



the companion because it neither manifests to us in the visible nor in the radio 
regions of electromagnetic spectrum. We will assume for it 10~^ s, which is rather 
reasonable since the periods of the known pulsars range from 1.56 x 10"'^ s to 
6 X 10"^ s. An-order of-magnitude calculation^ with our formulas yields 

/^\ -(4.6 X 10"'^ deg/y) X sin 2u;. (31) 

\ at / pgj^ B1913+16 

^Notice also that from pulsar orbit data reductions the inclination i cannot be determined unam- 
biguously. 

^It is mainly due to the classical effects of the oblateness of SS 2883; the relativistic gravitoelectric 
Einstein's precession amounts to only 3x10"^ deg/y. 

^The spin of PSR B1913-I-16 may be affected by the relativistic geodetic precession^, but since 
its period would amount to 360 years it can be considered constant over an orbital revolution. 
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Notice that the period of the secular precession of the perigee of PSR B1913+16 
amounts to 90 years^, so that over reasonable observational time spans of few 
years the predicted harmonics would be similar to secular trends. Since for PSR 
B1913+16 the experimental accuracy in measuring (w) amounts to about 1 x 10"^ 
deg/y (see reference^"), in this case we are not too far from the the possibility of 
detecting the predicted effect. 



5. Conclusions 

In this paper, in the linearized approximation of gravitoelectromagnetism and at 
order 0{c~^), we have calculated the influence of the spin s of a particle on its 
geodesic orbital motion in an external gravitomagnetic field generated by a rotat- 
ing body with spin J. It has been assumed that the orbit lies in the equatorial plane 
of the central object and that the spins are both vertical to it. The orbital-averaged, 
long-term effects on the Keplerian orbital elements of the orbiter have been calcu- 
lated by neglecting all terms of order C(e^). It turns out that the eccentricity, the 
pericenter and the mean anomaly are affected by such spin-spin gravitomagnetic in- 
teraction by means of harmonic perturbations with half the period of the pericenter. 
In order to see if they are measurable, we have examined two possible astronomical 
scenarios: the Sun-Mercury system and a pair of binary millisecond pulsar systems. 
It turns out that they are far too small. Only for the pulsar PSR1913-I-16 the pre- 
dicted periastron rate is not too far from the present experimental sensitivity. The 
possible discovery of new, highly eccentric and close binary pulsar systems, together 
with notable improvements of the accuracy in measuring the periastron rates could 
give some hopes to detect such tiny effects in future. 
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